Quantum criticality in a double quantum-dot system 
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We discuss the realization of the quantum-critical non-Fermi liquid state, originally discovered 
within the two-impurity Kondo model, in double quantum-dot systems. Contrary to the common 
belief, the corresponding fixed point is robust against particle-hole and various other asymmetries, 
and is only unstable to charge transfer between the two dots. We propose an experimental set-up 
where such charge transfer processes are suppressed, allowing a controlled approach to the quantum 
critical state. We also discuss transport and scaling properties in the vicinity of the critical point. 
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Quantum dots can be used to build single-electron 
transistors [1] and spin-based quantum bits [2], but 
equally interestingly, they serve as artificial atoms and 
allow to access correlated states of matter [3-5]. So 
far, most experiments focused on the study of Fermi- 
liquid states, with regular thermodynamic and transport 
properties at low temperatures [3, 5] and simple transi- 
tions or crossovers between them [4]. However, artificial 
molecules and mesoscopic structures can be used to real- 
ize and study non-Fermi liquids as well, characterized by 
singular properties and providing the simplest examples 
of quantum critical systems. However, due to their sin- 
gular nature, these states are very elusive. In fact, only 
recently Oreg and Goldhaber-Gordon [6] proposed a con- 
trolled set-up to access the two-channel Kondo (2CK) 
fixed point [7, 8], being the paradigmatic example of 
non-Fermi liquid impurity system. Subsequently, this 
setup was successfully realized experimentally [9]. Dis- 
sipation has also been proposed to drive quantum phase 
transitions (QPT) in quantum dots [10, 11]. However, 
most dissipative QPT are of Kosterlitz-Thoulcss type, 
and therefore no true quantum-critical state is realized. 

A non-Fermi-liquid state, similar to the one of the 
2CK model, emerges in the two-impurity Kondo model 
(2IKM). This model, initially studied in the context of 
heavy-fermion QPT, consists of two impurity spins that 
are coupled to conduction electrons and, at the same 
time, interact with each other through an exchange in- 
teraction. Jones et al. [12] observed that in the 2IKM a 
quantum critical point (QCP) separates a "local-singlet" 
from a Kondo-screened phase. This QCP has been shown 
to be essentially equivalent to the 2CK fixed point [14], 
though its operator content and finite-size spectrum are 
different [13]. In fact, it has been observed that - un- 
like the 2CK fixed point - the QCP of the 2IKM is 
very sensitive to certain electron-hole symmetry-breaking 
processes, which can smooth the QPT into a cross-over 
[13, 15]. (A related non- Fermi liquid fixed point also ap- 



peared in a two-orbital Anderson model [16].) 

The purpose of the present paper is to demonstrate 
that the QCP of the 2IKM can be realized and stud- 
ied in a system of two quantum dots, shown in Fig. 1. 
Such a double-dot system has a number of interesting 
regimes [17], however, here we shall focus on a situation 
far from the charge degeneracy points, with one unpaired 
electron on each of the dots. Remarkably, the quantum 
critical state in this geometry is very robust against both 
the asymmetry of the device (parity) and electron-hole 
asymmetry, and a sharp phase transition appears as long 
as there is no charge transfer between the dots 1 and 
2. We show that these charge transfer processes can be 
suppressed by inserting an artificial "antifcrromagnetic 
insulator" between the two dots (see Fig. lb). 

Model. To start our analysis, let us assume that the 
charging energies Eci,2 (Eci ~ Ec2 ~ Ec), associated 
with putting an extra electron to one of the two dots, 
are large compared to the level widths of the dots, T 7 
(7 = 1,2), and to the tunneling t between the two dots. 
Perturbatively integrating out virtual charge fluctuations 




FIG. 1: a) System of two quantum dots studied in the paper, 
b) Modified set-up with suppressed charge transfer processes, 
with an even number of quantum dots inserted between the 
two main dots (1,2) attached to leads. 
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of the dots, we arrive at the following Hamiltonian: 

= K 5i-5 2 + \{J$A + J^Si)^^, (1) 
+ {Vyji + Q 77 '5i ■ S^V^Vy + irrelevant terms. 

Here V 7<T = J2e 4yo-: witn c ! 7 tr being the creation 
operator of an electron state with spin a and energy 
e in the even combination of electrons in the leads at- 
tached to dot 7, and q 1 their density of states at the 
Fermi energy. Apart from irrelevant terms, Eq. (1) is the 
most general Hamiltonian that describes the double-dot 
system in the regime where charge fluctuations are sup- 
pressed. The largest couplings are K , J 7 = J 77 , and 
since these couplings are generated by second- 
order tunneling processes. They are typically of the size 
Ji ~ Vi ~ Ti/Ec, J 2 ~V 2 ~ T 2 /E c , and K ~ i 2 /£ c . 
The couplings Vi and V 2 can be made small by tuning 
the dots close to the middle of their respective Coulomb 
blockade valleys. The second-largest couplings are asso- 
ciated with charge transfer between leads 1 and 2, and 
are all of order V 12 ~ Qvx ~ J$ ~ [JiJ 2 K / E c ) 1/2 ■ 
All other couplings are suppressed by further powers of 
t/Ec, T/Ec, and do not change the physics essentially. 

Let us first study the Hamiltonian with the leading 
terms only, and no charge transfer between the two sides: 

flint = K S 1 S 2 + \{JiS 1 <4>{a^ l + JaSWIW (2) 

This Hamiltonian is characterized by three energy scales: 
Without the coupling K, the two spins on the two dots 
are screened independently at the Kondo temperatures 
Ti w Se e" 1 /-' 1 and T 2 « Se e' 1 ^ 2 , with Se < E c the 
typical level spacing on the dots [18]. These Kondo scales 
compete with K that tends to bind the two spins into an 
inter- impurity singlet. 

Clearly, the terms in Eq. (2) may break both parity and 
electron-hole symmetry. Nevertheless, solving Eq. (2) us- 
ing a numerical renormalization group (NRG) approach 
we find a sharp QPT upon variation of K for any value 
of the couplings J 7 and V 1 (in contrast to earlier state- 
ments). In all cases, the spectrum at the critical point 
can be described through a generalized version of the 
conformal field theory (CFT) of Affleck et al. [13], to be 
discussed below. 

Asymmetric limit. Before diving into the CFT solu- 
tion, let us give a simple and revealing physical picture 
of the physics in the limit T\ ^> K ^> T 2 . Here, the first 
spin is screened at a temperature T ~ T\. Below that 
scale, a local Fermi-liquid description applies to the re- 
sulting Kondo-screened complex, and therefore, it acts as 
a bath which tries to screen the spin S 2 [19]. The effective 
dimcnsionlcss coupling between S 2 and the Kondo com- 
plex can be estimated as Ai ~ KjT\. However, S 2 also 
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FIG. 2: Left: Phase diagram of the double-dot device in 
the absence of charge transfer. The two phases are sepa- 
rated by line of second order QPT, being very similar to the 
two-channel Kondo state. Right: Sketch of the temperature- 
dependent conductance through dot " 1" for Ti ~ T% ~ Tk , 
in the absence of charge transfer between the two sides. 

couples to spin excitations in the leads attached to it, 
with a rcnormalized coupling A2 ~ l/ln(Ti/T2). Clearly, 
we end up with an effective 2CK model, which is known 
to display a QPT at Ai = A2, corresponding to the con- 
dition T 2 w Ti cxp(— a Ti/K), with a a constant of the 
order of unity. The above argument is independent of 
the potential scattering terms. It shows that (i) The 
quantum-critical state is essentially identical to the two- 
channel Kondo state; (ii) Particle-hole or device (parity) 
symmetry are not required; (iii) The critical point is de- 
stroyed once there is charge transfer between channels 1 
and 2. The phase diagram obtained from these simple ar- 
guments is shown in Fig. 2. A similar picture is obtained 
within a CFT approach [13]. 

Conformal field theory. Since we do not have electron- 
hole symmetry in any of the channels, we used only the 
symmetries Ui(l) and ^2(1) associated with charge con- 
servation in the two channels and the global spin SU (2) 2 
symmetry for the conformal field theory solution. In 
the corresponding coset construction, Ui(l) x U 2 (l) x 
SU(2) 2 x Z 2 [13], all primary states and primary fields 
are characterized by their two charges, Qi an( i Q2, their 
spin j, and an Ising quantum number q (Id, a, and e). 

At the critical point the entire finite-size spectrum can 
be characterized just by two phase shifts, di,5 2 £ [0, 7r/2] 
. Similar to Ref. 13, the finite size spectrum is obtained 
through fusion with the Ising field a. The leading rel- 
evant operators at the fixed point are listed in Table I, 
where we also indicated the total charge, Q = Q\ + Q 2 
of every operator. Only operators with Q = can occur 
in the Hamiltonian, and in the absence of magnetic field 
only spinless operators can appear, therefore there are 
only two possible relevant operators, </> and <fr± that can 
be present in the Hamiltonian. Therefore, in the vicinity 
of the QCP, the Hamiltonian can be written as 

H = H* + + 5 c/) + + 5* (j) , (3) 
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TABLE I: Operator content of the critical point. 



where H* denotes the fixed-point Hamiltonian. The co- 
efficient k as 8K j \JTk = {K — Kc)/\/Tk measures the 
distance to the critical point, with Tk ~ min{Ti, T2} be- 
ing the Kondo scale associated with the formation of the 
non-Fermi liquid state. From the quantum numbers it is 
clear that the operators cj>± transfer exactly one charge 
from one side to the other, therefore the coefficient of 8 
is related to the amplitude of those operators in Eq. (1) 
that transfer charge between the two sides, and which 
have been neglected in Eq. (2). Both operators have scal- 
ing dimension 1/2 [20], and are thus relevant at the fixed 
point. However, k can be tuned to zero, while S always 
takes a finite value and generates a smooth cross-over at 
an energy scale Tj- ~ |(5| 2 , even for k = 0. As a result, a 
double-dot system never displays a true QPT. Neverthe- 
less, as we shall see later, the parameter 6 can be made 
small in a controlled way, such that the structure of the 
quantum critical point k = 8 = can be explored. 

Renormalization group. To obtain an estimate for the 
(dangerous) coupling 6 in Eq. (3) we need to compute the 
renormalization of the various processes that correspond 
to charge transfer in Eq. (1) [21]. To this purpose, let 
us assume that T\ as T2 as Tk and construct the pertur- 
bative scaling equations for the couplings in Eq. (1). In 
leading logarithmic order they read 



dl 



dQ 



(7) 



dl 



dV^ 
dl 
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Here I = ln(<5e/A) denotes the logarithmic energy scale, 
and we introduced a matrix notation in the lead indices, 
Q 7 y — > Q, . . .. From these equations we readily see that 
the most dangerous operators are the off-diagonal parts 
of the J™ which increase along the RG flow. However, 
in the perturbative regime the ratios / J\ and / J 2 
remain approximately constant. At the scale Tk wc have 
Ji ~ J2 ~ lj from which we immediately obtain an es- 
timate for the parameter 8: 8 ~ \/Tk{K / Ec) 1 ' 2 ■ Thus, 
for a double-dot system we find: Tg DD ~ TkK/Ec ■ For 
typical semiconductor quantum-dot parameters, Eq ~ 
20 K, and K ~ Tk ~ 0.5 K, this gives a cross-over scale 
T£ DD ~ 12 mK, which, while not very large, might be 
enough to spoil an observation of the non-Fermi liquid 
behavior. 



Suppressing charge transfer. T% can be suppressed by 
creating an artificial antifcrromagnctic insulator to medi- 
ate the exchange interaction between the two main dots 
1,2. The simplest arrangement is shown in Fig. lb, where 
we connect the two dots with two additional quantum 
dots with one electron on each of them. For simplicity, 
let us assume that the charging energies of all dots are 
similar, but the tunneling-generated exchange coupling 
K2 between the two central dots is somewhat larger than 
the one between the outer dots and their neighbors, Eq > 
K2 > K\ (see Fig. 1). In this limit, at energy scales below 
K2 the spins on the central dots are bound to a singlet, 
and their role is essentially restricted to mediate an anti- 
ferromagnetic interaction K ~ K 2 /K 2 between the two 
main dots. With parameters K 2 « 3 K and Ki as 1.5 K 
this gives a coupling in the range of K ~ IK ~ Tk- On 
the other hand, J$ ~ (JiJ 2 K 2 K 2 /E c ) 1/2 , and there- 
fore T$ is reduced to 



T* 



~ Tk[ e^) e- c 



(5) 



With the above parameters we find T$ 4D ai 10~ 3 Tk ~ 
0.5 mK. This value can readily be decreased even further 
by inserting more quantum dots in the middle. 

Transport. In the remainder of the paper we thus 
assume that T£ is smaller than the experimentally rel- 
evant temperature scales, i.e., we set 8 = 0. Let us 
furthermore concentrate on T\ as T2 ~ Tk- CFT al- 
lows to predict various observables in the regime close 
to the QCP, k as 0. We first note that in the absence 
of charge transfer, the linear conductance through dot 
7 is simply related to the T-matrix T^ of the conduc- 
tion electrons in the corresponding electrodes as G\ = 

($> im{T«/2} with g« = 4r xl r m /(r L1 + r R1 ) 2 

(see Fig. 1). At the fixed point (i.e., zero temperature), 
TW = i(l - SW), with S M the S-matrix of the elec- 
trons in lead 7 [22]. Similar to the analysis of [22] 
we find that = = at the QCP, and thus 

the conductance is G\(T = 0) = G^V 2 for K = K c- 
The approach to this value is determined by the leading 
irrelevant operator, which, similar to the electron-hole 
symmetrical case, is 0', the derived field from <j> [23]. 
At K = Kc, the finite-temperature corrections to Gq 
can be computed by perturbation theory in </)', with the 



result G hQC p (T) = G y Q L> [1 - a x ^T/Tk + ■■■)■ Here 
<X\ is a non-universal constant of order unity that de- 
pends on the asymmetry of the device and on the phase 
shifts. At finite source-drain voltages, V, the deviation 



8G\ = Gq — G±(T) will display scaling properties, sim- 
ilar to those of the 2CK model [24, 25] 



SGi/Gq = VT/T^ F{V/T) , 



(6) 



where the (non-universal) function F has the properties 
F(x <l)w const and F(x ^> 1) oc \fx. 
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FIG. 3: NRG results for Gi(w) (in unit of Go) for different 
RKKY couplings. The parameters used here correspond to 
Kondo couplings Ji ~ 0.1, J2 ~ 0.2, the potential scattering 
terms Vi ~ 0.003, V2 ~ 0.02, where the energy unit is the half 
bandwidth of the conduction electrons. The critical RKKY 
coupling is K c « 0.000976. The frequency u> in the plot is in 
units of Tk where Tk is defined as the half- width of Gi(w) 
for K = K C . 

For small but finite K, another crossover occurs at an 
energy scale T* = k 2 w (K - K c f /T K : For k > a 
inter-impurity singlet state is formed, while for k < a 
Kondo state is recovered. At these fixed points the S- 
matrices are given by S^ 7 ' = e 2%5 ~t (K > K c ) and = 
_ e 2i s., ^ < x c ), with both of these fixed points are of 
Fermi-liquid type, and therefore the conductance at them 
scales as Gi, sing iet = G (sin 2 (<5i) + p 1 (T/T*) 2 + . . .) and 
Gi.scrocncd = G (cos 2 (5i) - 71 (T/T*) 2 + ...), respec- 
tively, with fli and 71 again non-universal constants of 
order of unity [26] . The properties of G\ (T) are summa- 
rized in Fig. 2. 

A numerical computation of the finite-temperature 
scaling functions in the vicinity of the QCP is notoriously 
difficult. However, we can compute the AC conductance 
Gi(w) [27] by applying the NRG approach to the Ander- 
son Hamiltonian corresponding to Eq. (2). The results of 
this calculation for a generic situation without particle- 
hole and parity symmetries are shown in Fig. 3. The 
various crossovers can be clearly observed in Gi(uj) as a 
function of frequency, which displays a behavior qualita- 
tively similar to Gi(T). 

Summary. We have demonstrated that the quantum 
phase transition of the two-impurity Kondo model can 
be experimentally accessed using double quantum-dot 
devices. The non-Fermi liquid state is robust against 
particle-hole and device asymmetries; it is destroyed by 
charge transfer between the two main dots, which, how- 
ever, can be effectively suppressed with additional quan- 
tum dots in the set-up. Using a combination of analytical 
and numerical methods we have made predictions for rel- 
evant energy scales and transport quantities. 
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